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ABSTRACT 

In this paper we examine the contribution of galaxies with different infrared (IR) 
spectral energy distributions (SEDs) to the comoving infrared luminosity density, a 
proxy for the comoving star formation rate (SFR) density. We characterise galaxies as 
having either a cold or hot IR SED depending upon whether the rest-frame wavelength 
of their peak IR energy output is above or below 90 /j,m. Our work is based on a 
far-IR selected sample both in the local Universe and at high redshift, the former 
consisting of IRAS 60 /im-selected galaxies at z < 0.07 and the latter of Spitzer 70 /im 
selected galaxies across 0.1 < z ^ 1. We find that the total IR luminosity densities for 
each redshift/luminosity bin agree well with results derived from other deep mid/far- 
IR surveys. At z < 0.07 we observe the previously known results: that moderate 
luminosity galaxies (Lj_r < 10 11 Lq) dominate the total luminosity density and that 
the fraction of cold galaxies decreases with increasing luminosity, becoming negligible 
at the highest luminosities. Conversely, above z — 0.1 we find that luminous IR galaxies 
(Lir > 10 11 Lq), the majority of which are cold, dominate the IR luminosity density. 
We therefore infer that cold galaxies dominate the IR luminosity density across the 
whole < z < 1 range, hence appear to be the main driver behind the increase 
in SFR density up to z ~ 1 whereas local luminous galaxies are not, on the whole, 
representative of the high redshift population. 
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1 INTRODUCTION 

The rise in the comoving star formation ra te (SFR) density 
up to z ~ 1 and its subsequent flattening (|Lillv et al.lll996l : 
iMadau et all Il996l) has now been well studied at several 
wave lengths (e.g. iBunker et alj 120041 : iHopkins fc Beacoml 
120061 . and references therein). While the global picture to 
z = 1 has been well constrained by observation, the details 
of how or why this change occurs remain largely unknown. 
There now is evidence that star formation depends on both 
galaxy mass (iFeulner et alj 120051; lJuneau et al.l 120051) and 
environment l|Lewis et alj 20021 : Elbaz et al. 20071 ) present- 
ing a more complicated picture of galaxy evolution than sim- 
ple evolution of the luminosity function. 

The infrared (IR), and particularly the far-IR, is one 
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of the most powerful tracers of star formation a s IR lumi- 
nosity directly scales with SFR (|Kennicuttll 19981 . and refer- 
ences therein) and far-IR emission originates from regions 
of cold dust and gas that constitute the fuel for an on- 
going burst of star formation. Another advantage of select- 
ing sources at long wavelengths is the low contribution by 
active galactic nuclei (AGN), if present, to the total I R lumi- 
nosity l| Alexander et alj 12003 ; IClements et alj|2008l ). Stud- 
ies of the distant Universe at IR wavelengths remain lim- 
ited and the most sensitive probe has been surveys with 
the Spitzer 24 /im band. However, this wavelength is rela- 
tively far from the peak of all but the hottest IR galax- 
ies and progressively shifts to shorter wavelengths at higher 
redshifts where strong spectral features in the observed 
frame can also complicate matters. The few studies done 
with deep Spitzer 70 fim imaging have found rapid evolu- 
tion in the total IR luminosity function (|Huvnh et al.|[2007l ; 
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iMagnelli et"afl 12009). When the IR luminosity function is 
integrated at different redshifts it is found that the IR lu- 
minosity density increases rapidly up to z — 1 in a similar 
fashion to the SFR density. Nevertheless, the relative contri- 
bution by luminosity to the IR luminosity density changes 
with redshift with starbursts (10 ^ \og(Lm/ Lq) < 11) 
dominating locally, but with luminous IR galaxies (LIRGs: 
11 ^ \og(Lm/ Lq) < 12) and ultra-luminous IR galax- 
ies (ULIRGs: 12 ^ log(LiR/ Lq) < 13) becoming increas- 
ing contributors at higher redshifts (|Le Floc'h et al.ll2005l ; 
IMagnelli etalll2009l). 

In ISvmeonidis et aL (2009, hereafter S09) we studied 
the IR spectral energy distributions (SEDs) of a sample 
of 70 fim selected galaxies at z ^ 0.1, using a sample of 
local, z < 0.1, IRAS 60 /im selected galaxies for compari- 
son. We fitted the mid to far-IR photometry o f both sam- 
ples with models from ISiebenmorgen fc Kriigell (|2007T l. We 
found that the majority of the 70 fim sources had IR SEDs 
which peaked, in v x F v , at longer wavelengths than galaxies 
from the local sample with similar luminosities. In the local 
sample we observed a shift of the IR SED peak to shorter 
wavelengths with increasing luminosity (as has previously 
been noted: 
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witn increasing luminosity (as nas previously 
Sanders fc Mirabelll 19961 ; IChapman etaill2003l ; 



20091 ). whereas the 70 fim sample had a wide 



range of peak wavelengths the distribution of which varied 
little with luminosity. The observation that the IR SEDs 
of luminous, distant galaxies were on average different to 
their local anal ogue s was in contrast to other recent results, 
IMagnelli etail (e.g. 120091 ) who concluded there was no sig- 
nificant change in IR SED of luminous galaxies with red- 
shift. Models of galaxy evolution in the IR have tended to 
assume that the SEDs of high redshift lumin ous sources fol- 
low the luminosity trend seen in local sources llLagache et~al 
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2009). However, as shown in S09, the range of IR SEDs for 



luminous galaxies is much wider at high redshifts than seen 
locally. 

In this paper we examine the contribution of galaxies 
with different IR SEDs to the comoving IR luminosity den- 
sity (IRLD). We will compare our results to earlier work 
|Le Floc'h et al.ll2005l ) by examining the contribution to the 
IRLD by luminosity, but this time with a sample selected at 
70 /im rather than 24 fim. This selection enables us to use 
a more robust selection of IR luminous galaxies as 70 fim 
lies closer to the peak of typical IR SEDs. Hence, the to- 
tal IR luminosity can also be estimated more accurately, 
especially with constraints (mainly detections) from 160 fim 
data. We examine the IRLD in more detail by focusing on 
the contribution by IR SED type within each bin. By ob- 
taining an estimate of the peak of the IR SED in rest-frame 
v x F„, we can characterise these IR bright sources as be- 
ing either cold or hot depending on whether the SED peaks 
above or below 90 fim (32.2 K for a blackbody). We choose 
this wavelength for two reasons. Firstly, local IR galaxies 
appear to have a warm IR component (A ~ 60 (im) asso- 
ciated with dust around young star forming regions and a 
cooler "cirrus" component (A ^ 100 (im) associated with 
more extended dust heated by the interstellar radiation 
field l|Lonsdale fc Heloulll987T l. Secondly, this 90 fim cut also 
marks a divide between the SEDs of most local ULIRGs 
and those of cold SMGs seen at high redshif t, most of which 
would be classified as cold by our definition l|Chapman et al.l 




log(l + z) 

Figure 1. The redshift/total (8 — 1000 fim) IR luminosity dis- 
tribution of the local (z < 0.07) BGS sample and our Spitzer 
70 fim sample (z > 0.1) with sources characterised as cold or hot 
depending on whether their IR SED peaks above or below 90 fim 
(rest-frame in v X F v ). The dotted lines indicate the binning cho- 
sen to investigate the comoving IRLD by IR luminosity and SED. 



2003!) . We note that there is a linear relationship between the 
IRAS colour s often used to define the IR galaxies as 'cold' 
or 'hot' fe.g. lChapman et al.ll2003l ) and the peak wavelength 
used in this work. 

We present our far-IR galaxy sample in §2, describe our 
determination of the comoving IRLD in §3. We present our 
results in §4 and discuss them in §5. Throughout we use 
the current 'concordance' cosmology: Qm = 1 — = -3, 
Q = 1, and H = 70 krns" 1 Mpc" 1 (|Spergel et al.ll2003l ). 



2 THE FAR-IR GALAXY SAMPLE 

The IR luminous galaxies we use here are taken from the 
study of S09 who used a local and distant far-IR selected 
sampl e. The local sample i s the IRAS Bright Galaxy Sample 
from ISanders et al.l (|2003l ) which is selected to have 60 
flux densities above 5.24 Jy covering the entire sky surveyed 
by IRAS at Galactic latitudes |6| > 5°. The higher red- 
shift sample we use consists of Spitzer 70 fim galaxies from 
two well studied extragalact ic fields: the Ext ended Groth 
Strip (EGS covering - 0.3°. lDavis et alj|2007f ) and the 13 H 
XMM-Newton/Chan dra Deep Field (13 H covering ~ 0.3°, 
ISevmour et~ai1l2008l ). 

We restrict our 70 fim sample to sources with 0.1 ^ z ^ 
1.0 ruling out local resolved sources at z < 0.1 and very high 
redshift sources where there may not be complete redshift 
information. Although the 13 H field 70 fim sample has com- 
plete redshift information (40 spectroscopic and 40 photo- 
metric) the original EGS spe ctroscopic survey o nly targetted 
~ 2/3 of candidate sources llDavis et al. 2007f). By examin - 
ing the original EGS 70 fim sample ( Svmeonidis et ai]|2007l ) 
we estimate from broad-band colours that 101 sources lie 
within 0.1 ^ z ^ 1.0, 66 of which have spectra. Hence, in 
the next section we multiply the area of the EGS by a fac- 
tor 66/101 to correct for those sources not included on the 
assumption these 66 are representative of the parent pop- 
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Figure 2. The comoving IR luminosity density (IRLD) plotted as a function of redshift in different luminosity bins: starbursts (blue), 
LIRGs (yellow) and ULIRGs (red). For each bin we plot the total IRLD and the separate contribution from the cold and hot galaxies. 
The hot and cold points are marginally offset in redshift from the total IRLD for clarity. Derivation of the uncertainties and upper limits 
(representing null detections) is explained in the text. Locally (z<0.1) we use the IRAS Bright Galaxy Sample and at other redshifts 
we use the Spitzer 70 fim selected sample from S09. We overlay, for comparison, the IRLD derived from deep 24 fim observations by Le 
Floc'h et al. (2005, shaded regions). We observe a general trend of hot galaxies dominating the local L1RG and ULIRG bins, but the 
increase in IRLD at higher redshifts is mainly due to an increase in the contribution from cold galaxies rather than hot galaxies. 



ulation. We impose a further cut of <SVo/im > 10 mjy, to 
remove lower signal-to-noise sources where we are not com- 
plete, leaving us with 123 sources in our 70 fim sample from 
both surveys. 

We use the total IR luminosities and peak wavelengths 
dete rmined in S09 which were deriv ed from fits to models 
from lSiebenmorgen fc Kr iigcl (2007]). These models were fit 
to photometry at 12, 25, 60 and 100 fim for the local sample 
and at 8, 24, 70 and 160 fim for the Spitzer 70 fim sample. 
The IR luminosity-redshift distribution for both samples is 
plotted in Fig. [1] where we have characterised each source 
as having a cold or hot IR SED depending whether it peaks 
above or below our fiducial wavelength of 90 fim. 



3 METHOD 

We determine the IR lumi nosity density using a versi on of 
the classic 1/V max method (|Schmidtlll968l ; lFeltenlll976l ). We 
use the following formula to determine the luminosity den- 
sity in each redshift/luminosity bin for the 70 fim sample 

*m(£,«) = {Li/Vi) (1) 

where Li is the total 8 — 1000 fim IR luminosity (in solar 
units) of a source and Vi is the maximum volume (in Mpc 3 ) 
a source could occupy given the size of the redshift bins, the 
area of the survey, our lOmJy completeness limit and the 
best fit SED as determined in S09. The same method was 
used for the local BGS sample, also using the best fit SED 
found in S09, but the Vi was limited by the 5.24 Jy detection 



limit at 60 fim and a redshift range of 0.005 ^ z 0.07. 
Uncertainties are determined using the expression: 

8*m(L,z) = y/Ei (L^/V) 2 + 5®(8L l ) 2 (2) 

where Vi and Li are defined as before and the SQ(SLi) term 
represents the contribution of the uncertainty in the lumi- 
nosity of individual sources within that bin (due to photo- 
metric redshifts and the modeling in S09) . The binning (see 
Fig. [1} was defined by the definitions of the luminosity class 
(ULIRG etc.) and the desire to have bins of equal size in 
log(l + z) space, i.e. approximately equal epochs in cosmic 
time. We do not include a bin to cover the highest redshift 
LIRGs as our survey is much less sensitive to cold galaxies 
in this region of parameter space and hence we are unable 
to derive a representative ratio of hot and cold galaxies. For 
each luminosity-redshift bin we calculated the total comov- 
ing IR luminosity density as well as that for the cold and hot 
galaxies separately. The redshift used when plotting each bin 
is the mean redshift of all the sources contributing to that 
bin in log(l + z) space. Upper limits for a bin with no cold 
or hot sources were determined from the 97.7 2% confidence 
limit (equivalent to 2 a) for a null detection (lGehrelslll986l ) 
using the volume of the whole bin and the mean luminosity 
of all the galaxies in the bin. 

In S09 we discussed how the construction of our sample 
might affect the results. The selection of the 70 fim sample at 
rest-frame 35 — 60 fim means our selection could potentially 
favour galaxies that peak around those wavelengths, i.e. hot 
galaxies by our definition, contrary to the increasing fraction 
of cold galaxies. We also established that the uncertainty 
in the long- wavelength Spitzer/ 160 fim photometry did not 
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Figure 3. The fraction of the comoving IR luminosity density for 
different redshift/luminosity bins due to cold galaxies (A pca k > 
90 (im) also given as a percentage on the right axis. Note by this 
definition we have no cold galaxies in our local ULIRG sample. 



have a significant effect on the peak wavelength determined 
from the SED fitting. The uncertainties on IR luminosity 
from the SED fitting, when carried over to this work, con- 
tribute little to the total uncertainties on individual data 
points in Fig. [5] i.e. small number statistics dominate. 



4 RESULTS 

We present the observed comoving IR luminosity density of 
our sample in Fig. [5] as a function of luminosity and redshift. 
The total IR luminosity densities at each luminosity-redshift 
bin are in good agreement with those from Le Floc'h et al. 
(2005) notwithstanding the potential sample variance be- 
tween small volumes and the different IR SED models used. 
The results presented here also broadly agree with those of 
Magnelli et al. (2009) who base their work on stacking 24 /im 
selected sources at longer wavelengths. The largest disagree- 
ment with previous work is at z ~ 0.2 where the total IR 
luminosity density of the starbursts presented here is about 
0.4 dex below the value of Le Floc'h et al. (2005). 

Examining the contribution by IR SED to each 
luminosity-redshift bin, we find that there is significant vari- 
ation in the relative contribution of cold and hot galaxies 
with luminosity and redshift. Locally we see cold sources 
dominating the starburst bin, but at higher luminosities 
they contribute less; cold galaxies comprise only one third 
of the LIRG bin and they contribute nothing to the ULIRG 
bin. The contribution of cold galaxies to the IR luminosity 
density from starbursts remains fairly constant at 60 — 70% 
over the redshift range we can detect them, ^ z ^ 0.2. 
However, for both LIRG and ULIRG bins the contribution 
from cold galaxies increases considerably at higher redshifts. 
Cold galaxies become the dominant contributor to the IR lu- 
minosity density in these bins above z = 0.1, in contrast to 
the situation seen in the local Universe. In fact, cold galaxies 
are responsible for most of the rise with redshift of the IR 
luminosity density in both the LIRG and ULIRG bins. 

In Fig. [3] we plot the fraction of the total IR lumi- 
nosity density in each redshift/luminosity bin due to cold 



galaxies. We can see how the contribution of cold galaxies 
to the starburst luminosity density remains constant up to 
2 ~ 0.2, although our data here cannot constrain that frac- 
tion at higher redshifts. The contribution of cold galaxies to 
the LIRG and ULIRG luminosity density increases dramat- 
ically with redshift; locally the contribution is very low or 
negligible, but at higher redshift cold galaxies dominate. We 
find that our results do not change qualitatively if we change 
our choice of fiducial wavelength (that divides between hot 
or cold galaxies) by ±10 /im. 



5 DISCUSSION 

In S09 we demonstrated that the IR SEDs of individual 
LIRGs and ULIRGs at 0.1 ^ z 1.0 peak at longer wave- 
lengths, and span a wider range, than local galaxies at sim- 
ilar luminosities. Other changes in the properties of star 
forming galaxies have occured since z ~ 1. Star formation 
shifts from massive galaxies at high redshifts to progressively 
less massive galaxies at lower redshifts (e.g. ICowie et al.l 
1996: J uneau et al.ll2005l ). A shift is also seen in environ- 
ment of the most strongly star forming galaxies with star 
format ion tending to oc cur in over-dense environments at 
z ~ 1 l|Elbaz et al.ll2007l ). a reverse of the trend seen locally. 
While the dependence of star formation on stellar mass and 
environment inform us about the star formation history and 
potential or current galaxy mergers, information about the 
IR SED is a direct measure of the physical conditions in 
galaxies underlying major episodes of star formation. 

In this paper we have expanded the results of S09 by de- 
riving IR luminosity densities at different luminosities and 
epochs. We observe that the total IR luminosity densities 
derived from a 70 /im sample are broadly consistent with re- 
sults derived from shorter wavelengths (Fig [2j . The low IR 
luminosity density of the z ~ 0.2 starburst bin compared 
to other results is likely due to the effect of sample vari- 
ance over the small cosmic volume probed. We also separate 
our sample into galaxies that can be characterised as hav- 
ing cold or hot SEDs and study their contribution to each 
redshift/luminosity bin. We find two striking new results. 
Firstly, we observed a rapid change in the make up of LIRG 
and ULIRG bins: locally they are dominated by hot galax- 
ies, but at z > 0.2 they are dominated by cold galaxies. 
Secondly, cold galaxies appear to be the dominant contrib- 
utor to the total IR luminosity density over ^ z ^ 1.0; 
when starburst galaxies dominate at z < 0.2 we find that 
they are mainly cold and by the time LIRGs and ULIRGs 
begin to dominate at higher redshifts they are also mainly 
cold. Hence, a cold mode of star formation has dominated 
galaxy evolution since z ~ 1 and is likely responsible for the 
increase in the star formation rate density up to z = 1. 

While Fig. [2] does dramatically show the rise of cold 
galaxies in the LIRG and ULIRG bins at higher redshifts we 
note that hot galaxies are still found at those redshifts. In 
fact the results presented here do not rule out a rise with red- 
shift in the number (and luminosity density) of hot galaxies. 
Therefore, the distribution of types of SEDs remains quite 
wide at high redshift. Fig. [3] most strikingly demonstrates 
the evolution of the relative contribution of cold galaxies to 
each redshift/luminosity bin. The rise in the contribution of 
cold galaxies to the LIRG and ULIRG bins is most rapid be- 
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tween z ~ and 2 ~ 0.3. We now consider if AGN could be 
responsible for the change in the fraction of cold galaxies. As 
a powerful AGN would make the IR SED hotter the presence 
of an AGN can not explain the prevalence of cold galaxies at 
high redshift. AGN could, however, be the cause of the high 
fraction of hot galaxies locally. Studies have shown that the 
incidence of AGN are low in local LIRGs, 10 - 20%, (Pet- 
ric et al. , 2009, in press), bu t while it is higher in ULRIGs, 
~ 40% (|Farrah et all 120071 ) the 8 - 1000 IR luminos- 
ity remains dominated by star formation. We note that at 
higher redshift the incidence of AGN, and t heir contribution 
to total IR luminosity, als o remains low (|Alexander et all 
120031 ; IClements et aLlf2008l ). Hence, AGN are not a signifi- 
cant contributor to the IR luminosity of the sources studied 
here and therefore are not responsible for any of the trends 
with redshift. 

The prominence of cold galaxies discovered here might 
suggest that there could be a substantial contribution to the 
IR emission of luminous galaxies at high redshift from a cold 
cirrus component akin to that observed in local IR galaxies. 
Such an interpretation is in agreement with the observa- 
tion that many of t he distant, cold, SMG s are extended on 
scales of ~ lOkpc (|Chapman et al.l I2004T ) in contrast with 
the compact ULIRGs seen locally, and is also consistent with 
the change in the mean IR SED of luminous galaxies as in- 
ferred from the observed 70 fim to radio flux density rati o of 
high redshift star forming galaxies ( Se ymour et al.ir2009T ). In 
addition, changes in dust properties (such as opacity, grain 
distribution etc.) with redshift are also a possible cause of 
the larger spread in SED types at high redshift. 

Most galaxy evolution models (e.g. | Lagache~et al.l 
|Pearsonll2005l ; iLe Borgne et al.ll2009l ; I Rowan- Robinson 
assume that the IR SED is only dependent on IR luminos- 
ity, i.e. they typically use certain templates for star forming 
galaxies of a given luminosity, often using local galaxies like 
Arp 220 (which can be characterised as having a hot IR SED 
by our definition) for the most luminous population. This 
choice was often necessary as there were few constraints on 
any change with redshift until now. Models based on cooler 
SEDs could , for example, reconcile the flat or decreasing 
star formation rate density above z — 1 and the far-IR and 
sub-millimetre source counts. 

In conclusion, we have found that the star formation 
history of the Universe is dominated by cold galaxies across 
< z < 1 and that the hot luminous IR galaxies we see 
locally are not representative of the luminous galaxies that 
make up the bulk of the star formation at 0.1 < z < 1. 
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